Implementation of Spin Hamiltonians in Optical Lattices 
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We propose an optical lattice setup to investigate spin chains and ladders. Electric and mag- 
netic fields allow us to vary at will the coupling constants, producing a variety of quantum phases 
including the Haldane phase, critical phases, quantum dimers etc. Numerical simulations are pre- 
sented showing how ground states can be prepared adiabatically. We also propose ways to measure 
a number of observables, like energy gap, staggered magnetization, end-chain spins effects, spin 
correlations and the string order parameter. 



In Condensed Matter Physics, there are strongly cor- 
related systems of spins and electrons whose study is 
extremely difficult both analytically and numerically. 
These systems are of great practical interest since they 
may be relevant in some important instances like high- 
Tc superconductivity, to quote just one open problem of 
great relevance. These open problems have triggered a 
great number of models formulated by means of quantum 
many-body Hamiltonians like Heisenberg, t-J, Hubbard 
etc. Their quantum phase diagrams remain unknown for 
generic values of coupling constants, electron concentra- 
tion (doping) and temperature, although a great knowl- 
edge can be obtained in particular integrable models in 
one dimension. 

Two very important examples of these systems are 
quantum spin chains and ladders. Since the seminal work 
of Haldane P] quantum spin chains have been exten- 
sively studied as one of the simplest but most emblem- 
atic quantum many-body systems. According to Hal- 
dane, the one-dimensional integer-spin Heisenberg anti- 
ferromagnets have a unique disordered ground state with 
unbroken rotational symmetry and with a finite excita- 
tion gap in the spectrum, while half-integer antiferromag- 
nets are gapless and critical. This quantum many-body 
phenomena is different from the usual source of gaps in 
magnets, namely, single-ion anisotropy, which does not 
involve quantum correlation effects. Several theoretical 
developments jj, have helped to clarify the situation 
and there is now strong numerical evidence in support of 
Haldane's claim. 

Here we shall address the issue of implementing spin 
chains and ladders in an optical lattice. For concreteness, 
we focus on spin chains first. For integer spin, s = 1, 
there is a quantum Hamiltonian that contains all the 
relevant information pertaining the Haldane phase and 
exhibiting a rich phase structure. This is the so-called 
Quadratic-Biquadratic Hamiltonian (QBH) given by 
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Here, Si are spins s = 1 at lattice site i, /? is a relative cou- 
pling constant that parametrizes a family of local Hamil- 
tonians, and the sign of a determines the ferro or antifer- 



romagnetic regimes. The properties of the ground state 
without magnetic field, B — 0, are entirely determined by 
an angle, 9, such that a = \a\ cos{9) and a/3 = — |a| sin(6') 
[See Fig.^a)]. For 9 £ [— f , f ], the ground state belongs 
to the Haldane phase, with 9 = being the Heisenberg 
point and 9 = arctan(i) the AKLT point 2], which is of 
particular importance because it can be described with 
an exact valence-bond wavefunction. There are two crit- 
ical points on which the standard correlation length 
diverges: 9 = j is the Uimin-Lai-Sutherland (ULS) crit- 
ical point Q , where a phase transition occurs into a gap- 
less phase for 9 £ [f , §]; and 9 = —j corresponds to the 
Takhtajan-Babujian (TB) critical point Q displaying a 
second order phase transition into a dimerized phase for 
9 S [^j^, —j], thus gapped. It has been conjectured that 
a quantum nondimerized nematic phase also exists |5|| in 
the ferromagnetic region. 

Likewise, ladders of spin s = ^ |^ exhibit a rich quan- 
tum phase structure depending on their number of legs 
and couplings: for even legs, ladders are gapped and can 
be in a Haldane phase similar to the integer spin chains, 
whereas for odd legs they are gapless. Even-legged lad- 
ders can also be in quantum dimer phases located on the 
rungs, provided the vertical couplings are strong enough. 
Theoretically, spin ladders are regarded as a route to ap- 
proach the more complicate physics of two-dimensional 
quantum spin systems, as we increase the number of legs. 
For instance, a two-leg ladder is gapped and upon hole 
doping can serve as a toy model for studying supercon- 
ducting correlations. 

Experimental attempts to implement the QB Hamilto- 
nian in real materials date back to the sixties. The first 
experimental evidence of a biquadratic term was found in 
Mn-doped MgO Q and for 5 = 5/2 Mn++ ions in anti- 
ferromagnetic chains of MnO . The measured coupling 
constant was (3 = —0.05, too small and with a fixed sign. 
Thus, only the pure Heisenberg AF model and its im- 
mediate surroundings with a very small /? are of feasible 
practical implementation, and the AKLT point is rather 
far to be accessible. Experimentally, ladders can be re- 
alized by selective stoichiometric composition of cuprate 
planes in superconductors 0- However, residual non- 
vanishing inter-ladder couplings on the planes introduce 
disturbances which are difficult to control. 
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FIG. 1: (a) DifFerent phases of the ground state of Hamilto- 
nian Q, according to the "angle" 6 defined in the text, (b) 
Type of Hamiltonian Q as a function of the gradient of the 
electric field, A, for Uo = 0.75f/2. The solid line is obtained 
naturally, the dashed line is for the dual model obtained when 
working on the upper part of the spectrum, and the dotted 
line marks the location of the possible nematic phase. 




FIG. 2: (a) A ladder is the combination of two spin chains 
that interact with each other. Interactions along a chain and 
between legs can be different, (b) We can build a ladder 
combining a 3D lattice, l2{x), with an additional ID optical 
lattice, Ii{x), that has twice the period of the first one. This 
induces a tunable hopping J', different from the longitudinal 
one, J, and suppresses hopping between neighboring ladders. 



Finding an experimental setup for checking the vahdity 
and observing the several phases of the QB Hamiltonian 
and of spin ladders is considered a very important chal- 
lenge in the field. In this paper we propose to solve this 
joblem using cold atoms confined in an optical lattice 
As shown before, a Mott phase 0| of cold atoms in 
a lattice can be described using ferromagnetic spin s = ^ 
[TSini or s = 1 [HEl Hamiltonians. Here we describe 
how to access a wider family of models, including Hal- 
dane phases of antiferromagnetic s = 1 chains and s = ^ 
ladders. We also design a technique to prepare adiabat- 
ically the atoms in the ground state, an important task 
since these spins cannot be cooled. Finally we study how 
to detect the different spin phases, and to directly ob- 
serve correlation and excitation properties. This is in 
sharp contrast with standard experiments in condensed 
matter where one neither has a controllable implemen- 
tation of the QB Hamiltonian ^ nor a direct way to 
perform measurements that so far have been regarded as 
mere theoretical tools. 

Let us first consider how to engineer Hamiltonian 
using spin s = 1 bosons in an array of one-dimensional 
optical lattices. For a strong confinement and low den- 
sities, the effective Hamiltonian is the Bose-Hubbard 
model nil 

H = -J {a]^aic + al^aja) + J^^^^ + -Sj")a]„aja 
+ E ^ E i<%^3fa3s)H<%a,^a,^)- (2) 

S=0,2 j,Q,/3,7,5 

First of all, while the indices j and / run over the lattice 
sites, the Greek letters label the projection along the Z 
axis of either the spin of an atom (a, /3, 7, (5 = —1, 0, -1-1), 
or of a pair of them (cr — —2, —1, 0, 1, 2). Then the first 
term in the Hamiltonian is the single-particle hopping 
term and J is the tunneling amplitude to a neighbor- 
ing site. The second term models the interaction be- 
tween bosons in a single site: two bosons can only in- 
teract if their total spin is either or 2, because the 



state 5 = 1 is antisymmetric. Furthermore, the inter- 
action may be different for each value of the total spin. 
Both statements are summarized in the presence of spin- 
dependent interaction constants, Us, and in the tensors 



^^s ~ ('S'j ^k, 7; s, /3), which are the Clebsch-Gordan 
coefficients between the states \s = 1, 7) ® |s = 1, (3) and 
\S = 0,2; a). Finally, we have included effective elec- 
tric and magnetic fields, Ej and Bja, that can be engi- 
neered using Stark shifts and spatially dependent mag- 
netic fields, as in current experiments 

We will assume that the lattice has been loaded with 
one atom per site |l5j . and that the tunneling has been 
strongly suppressed, J <^ Us- With a perturbative cal- 
culation around states with unit occupation 0, 0, 
we obtain the QB Hamiltonian with constants 



J^Us 
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This result is valid only if the gradient of the magnetic 
field is small, — Bj\ ^ \Us\, and the gradient of 

the electric field is constant, A — — Ej, and not 

resonant with the interaction, \nUs ± A| ^ J, Vn e Z. 

In the absence of electric or magnetic fields the model 
reduces to that of 0| and ^ , and we are restricted to 
a fixed value of 0, typically in the ferromagnetic sector. 
However, with our tools it is possible to explore many 
other phases and achieve almost all values of 9 [Fig. 
The idea is to change the gradient of the electric field and 
use a duality H^p = —Hp between ferro and antiferro- 
magnetic models: The highest energy state of a ferromag- 
netic model (a < 0) is the same and exhibits the same 
dynamics as the ground state of the dual model {—a, —(3), 
since {idt- Hp)i}{t) = <^ (idt - Haf) V{t) = 0. 
This equivalence is possible in current experiments, be- 
cause dissipation is negligible and decoherence affects 
equally both ends of the spectrum. 

A similar procedure is used for implementing ladders 
of spin s = i . A ladder is nothing but the combination 
of two spin chains (legs) that interact with each other. To 
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FIG. 3: (a) Procedure for the adiabatic construction of the 
ground state of the s = 1 antiferromagnet with open bound- 
ary conditions, (b) Infidehty of the final state for the AKLT 
(dashed) , Heisenberg (solid) and s = | spin ladder (dots) , as 
a function of the duration of the process, for a setup with 9 
spins (chains) or rungs (ladders). 



build them we need to set up a 3D lattice that confines 
the atoms on planar square lattices (hopping has been 
suppressed along the Z direction) , and superimpose along 
the Y direction a second ID lattice with twice the period 
[Fig-lljb)]. Adjusting the intensities of different lattices 
we can modify the tunneling along the leg of a ladder 
and between neighboring legs, and completely suppress 
tunneling between ladders. With the help of electric and 
magnetic fields fT]l|, and the duality between ferro and 
antiferromagnetic models, we achieve once more a full 
tunability of the Hamiltonian. 

Let us now study how to prepare ground states adi- 
abatically. We will focus on the Haldane phase of the 
s = 1 lattice and on the antiferromagnetic s = ^ ladders. 
Since in both cases we seek an antiferromagnetic state, 
we can begin with a configuration of antiparallel spins, an 
effective staggered magnetic field, Bj = {—iy\B{t)\, and 
no hopping. We then progressively decrease the magnetic 
field and increase the interaction, a [See Fig. Ol^a)]. This 
procedure constrains us to a subspace of fixed magneti- 
zation, 5*1) ~ 0, and also ensures that the minimum 
energy gap between the ground state and the first excita- 
tions remains independent of the number of spins. Thus, 
the speed of the adiabatic process can be the same for all 
lattice sizes, an important point in a setup with defects. 

We have studied numerically the fidelity of the adia- 
batic process for the AKLT point, for /3 = 0, and for 
a s = ^ ladder, using different speeds and sizes. The 
fidelity is the projection of the final state onto the (de- 
generate) ground states (in the case of s = 1, it is one 
singlet and three triplet states). The results are shown in 
Fig. I^b) . Already for a duration of 10 x h/amax the Hal- 
dane phase is built with 99% fidelity. Assuming current 
optical lattices, with U /h ~ 3KHz and J = 0.3C/, this 
implies a time of roughly 37 ms. Improvements on these 
values are expected with the implementation of optically 
induced Feschbach resonances ,16.] . 

Regarding the staggered magnetic field, for ^^Rb in the 
-F = 1 hyperfine state it can be produced using a weaker 
optical lattice, aligned with the atomic chains and made 
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FIG. 4: Maximum staggered magnetization per spin ac- 
quired by the ground state of the (a) AKLT and (b) Heisen- 
berg models, under an oscillating magnetic field, Bj — 
(-1)^0. 025a cos(27ri/t), for 5, 7 and 9 spins (right to left). 
Vertical lines mark the finite excitation gap for the given size. 



of a pair of counterpropagating laser beams in a lin _L 
lin configuration. If the lasers are far off-resonance from 
the transition ^S ^P, we will obtain a state-dependent 
potential V±{x) oc 2 + (sin(fca;)^ — cos(A:a;)^)5'^, where 
the sin(A:a:) and cos(fca;) come from the Stark shifts in- 
duced by the (T+ and (t_ polarizations on the atomic 
states. Choosing the orientation of the counterpropa- 
gating beams so that V±{x) has twice the periodicity of 
the confining lattice we get our staggered magnetic field. 
A similar setup can be designed for s = i particles. 

Once we have constructed the ground state, we would 
like to study its properties. We will describe a number of 
possible experiments, sorted by increasing difficulty. For 
illustrative purposes, we consider the Haldane phase of a 
s = 1 spin chain, but we want to emphasize that the same 
techniques can be applied to other phases, spin models 
and even ladders. Preliminary evidences of the Haldane 
phase can be obtained by studying global properties, such 
as the staggered magnetization, Sst ~ J^ji^^V ^jy which 
is zero in the dimer phase and nonzero in the Haldane 
phase. To measure S^^^ , we apply a 7r/2 rotation around 
the Z axis using the staggered magnetic field, and then 
measure Q2iS^'^). The remaining component S^^ can 
be obtained by rotating all spins an angle 7r/2 around 
the X or Y axes, and then measuring S'^'j or Sf^.. 

We can also study the energy gap between the ground 
state and its excitations, using an oscillating magnetic 
field, Bj{t) = {—ly B sm{ojt) . From linear response the- 
ory we know that for small intensities, \B\ a, there is a 
strong resonance at the gap, huj = Egap, which manifests 
itself on the growth of the staggered magnetization, S'J^. 
This result has been confirmed by numerical simulations 
of small lattices, as shown in Fig. 01 

Another interesting feature of the Haldane phase are 
fractionalization effects. To understand this, one should 
visualize each atom with spin s = 1 as being composed 
of two s = ^ bosons in a symmetric state. The ground 
state of Eq. ^) can then be built — either approximately, 
if P ^ 1/3, or exactly, for the AKLT — , by antisym- 
metrizing pairs of virtual spins from neighboring sites 



P|. This leaves us with two free effective s = ^ spins at 



4 



the ends of a chain, which manifest themselves physically. 
First, the four almost degenerate ground states are de- 
termined by the values of the free virtual spins, which we 
will denote as |«i,?2Af)- Thus, if the state of the system 
is = Cij the probability of measuring the 

left- and rightmost real spins in states I, J = ±1 is ap- 
proximately 4|c//2,j/2p/9- And second, the virtual spins 
almost do not interact and can be manipulated indepen- 
dently with weak magnetic fields that have different val- 
ues on the borders of a chain. For instance, if we prepare 
the ground state using our method, apply a global rota- 
tion of angle ^ = f around the Y axis, and then switch 
on the staggered magnetic field, we will measure periodic 
oscillations in the value of Sff , because the virtual spins 
on even and odd sites rotate with opposite senses. 

We also have developed a procedure to measure spin 
correlation functions in optical lattices, an essential tool 
for experiments with spin lattices. Our proposal only as- 
sumes that we can trap atoms in separate lattices and 
empty lattice sites with double occupation |23|. First 
of all, we notice that a spin correlation may be written 



as a density correlation, (S^S^) 



+1 -1 -1 +: 



where is the number of atoms in 



hyperfine state a on the A:-th site. A correlation such 
as ('^}''^}+a) *^an- be measured by moving the atoms 
of species +1 just A sites 0, emptying all doubly 
occupated sites and counting the number of atoms left 
in state -f 1, on the sites j and k. If we rather count the 
atoms in state 0, then we will obtain {nj^n^^). With 
this method, and some global rotation of the spins, it is 
possible to obtain all correlators, (SjSj+A)- If we can- 
not address individual atoms, using the same procedure 
and measuring total populations, X^fc'^fej ^^^^ obtain 
averaged values, J2ki^j^j+^)- Both quantities are inter- 
esting to discriminate between the different spin phases. 

On a much higher level of difficulty, but with pretty 
much the same tools as for measurin g co rrelations we can 
obtain the string order parameters ly|. One can show 
that it is equivalent to a correlation function measured 
on a transformed state, 



5,™ := (^,^e-S7-+i^^5^,)^ = (5,^5;;,)^^_^^(4a) 



Ua — exp 



N k+A 
k=lj=k+l 



The unitary operation Ua can be performed as fol- 
lows. First we perform half a swap between the -1-1 and 
-1 states, C/i = exp [^^2 ^^(1+1) (-11^ + 1-1) J], 
with a 7r/2 Raman pulse that connects these states. Next 
we split the three species into three optical lattices. The 
atoms in state will move A — 1 sites to the right, and on 
each movement a controlled collision with atoms in state 
-1-1 will take place. Adjusting the duration of this colli- 
sion so that it produces a phase of tt, we obtain the trans- 



formation U2 = exp 



(4b) 



iEkEfJiim {o|J+i)(+iL.+,) 

We restore all atoms back to their positions, and repeat 
the operation [/i, concluding the total transformation 
Ua — U1U2U1. Finally, we perform all required steps 
to measure either (S^S'^+a) or T,k{SkSk+A)- 



Typical experiments have defects and thus host chains 
with different number of spins. However, except for the 
string-order parameter, the measurements that we pro- 
pose are extremely robust, and in general they produce a 
signal that is a nonzero average of the possible outcomes 
for chains of different lengths. 



Summing up, in this work we have shown how to im- 
plement spin s = 1 chains and s — ^ ladders with cold 
atoms in an optical lattice. Such experiments will allow 
us to construct never observed phases and probably dilu- 
cidate the existence of the nematic phase. Finally, we 
have developed a very general set of tools to characterize 
these spin phases, which are themselves of interest for 
future experiments with optical lattices. 
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